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ABSTRACT 


Improved  methods  are  presented  for  estimating  the  required  start¬ 
ing  pressure  ratio  of  a  given  axisymmetric  second-throat,  ejector- 
diffuser  system  used  in  rocket  altitude  simulation.  The  improved  theory 
involves  a  more  accurate  technique  for  determining  ramp  pressure.  The 
new  technique  utilizes  the  flow  field  approaching  the  second  throat  along 
with  two-dimensional  oblique  shock  relations  and  the  conservation  equa¬ 
tions  to  determine  the  conditions  at  the  throat  entrance.  The  results 
are  applied  to  second-throat  diffusers  with  throat  lengths  varying  from 
0.  2  to  8.  1  throat  diameters  and  to  second  throats  with  a  centerbody 
located  at  various  positions.  A  comparison  of  the  theoretical  starting 
pressure  ratios  with  the  experimental  data  from  previous  work,  as  well 
as  the  additional  experimental  data  obtained  during  this  investigation, 
gave  a  maximum  deviation  of  less  than  10  percent  for  most  of  the  con¬ 
figurations  investigated.  The  improved  method  should  be  especially 
applicable  to  ejector-diffuser  systems  whose  nozzle  exit  flow  fields  are 
nonuniform,  such  as  exist  in  annular  and  clustered  driving  nozzles. 
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The  application  of  second-throat  diffusers  in  testing  of  rocket  motors 
at  simulated  altitude  conditions  has  gained  an  increasing  interest  at  the 
Rocket  Test  Facility  (RTF).  With  the  increase  in  diffuser  pressure  recov¬ 
ery  obtained  from  a  second-throat  diffuser,  higher  altitudes  can  be  simu¬ 
lated  by  the  ejector-action  of  the  rocket  exhaust  gases  with  the  same 
facility  exhaust  machinery. 

The  investigation  of  second-throat,  ejector -diffuser  systems  in  RTF 
was  reported  initially  in  Ref.  1.  Semi-empirical  methods  were  later 
developed  in  Ref.  2  to  calculate  the  diffuser  starting  pressure  ratio  which 
gave  very  good  agreement  with  the  experimental  results  for  both  long 
[  (L/D)s^  >  5  ]  and  short  [  (L/D)st  <  1]  second-throat  ejector  systems 
using  conical  driving  nozzles.  The  short  second-throat  theory  utilizes 
the  conservation  equations  to  arrive  at  a  one -dimensional  supersonic 
Mach  number  at  the  entrance  of  the  minimum  area  at  station  2  (see  Fig.  1). 
The  maximum  back  pressure  was  then  determined  which  would  cause 
separation  of  a  turbulent  boundary  layer  based  on  the  free,  shock - 
separation  criteria  proposed  by  Arens  and  Spiegler  (Ref.  3).  The  long 
second-throat  theory  used  the  same  conservation  equations  to  arrive  at 
the  force  balance  between  stations  1  and  2.  However,  the  subsonic  solu¬ 
tion  was  used,  together  with  the  continuity  equation,  to  calculate  the  maxi¬ 
mum  diffuser  exit  pressure  at  station  3  for  which  the  long  second-throat 
diffuser  will  start.  Both  of  these  theories  involved  a  semi-empirical 
method  for  predicting  the  pressure  drag  force  on  the  second-throat  ramp, 
which  assumed  that  the  gas  expands  isentropically  from  the  nozzle  throat 
to  the  diffuser  duct  diameter  through  a  mean  gas  expansion  angle,  0m,  as 
shown  in  the  following  sketch. 


This  procedure  provides  a  simple  and  sufficiently  accurate  method  for 
predicting  the  conditions  on  the  ramp  and  approaching  the  second-throat 
entrance  for  second-throat  ejector  systems  using  conical  nozzles  or  noz¬ 
zles  with  uniform  exit  Mach  numbers.  However,  this  semi-empirical 
method  is  not  adequate  for  annular  or  cluster  nozzles,  whose  flow  condi¬ 
tions  at  the  nozzle  exit  are  highly  nonuniform.  The  investigation  was 
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therefore  continued  to  find  a  more  general  approach  for  calculating  the 
starting  pressure  ratio  of  a  given  second-throat  ejector  system.  A 
satisfactory  theoretical  method  for  determining  the  performance  of 
such  systems  was  developed,  which  requires  that  the  flow  conditions 
approaching  the  second  throat  be  known.  The  pressure  drag  force  on 
the  second-throat  ramp  and  the  conditions  at  the  throat  entrance  are 
then  determined  by  using  two-dimensional  oblique  shock  relations  and 
the  conservation  equations.  The  starting  pressure  ratio  for  a  given 
second-throat  diffuser  configuration  can  then  be  calculated  by  making 
additional  assumptions  for  the  forces  imposed  downstream  of  the  throat 
entrance,  depending  on  the  type  of  second-throat  diffuser  under 
investigation. 

Before  discussing  second-throat  diffusers,  it  is  necessary  to  define 
what  is  meant  by  starting,  operating,  and  minimum  cell  pressure  ratios, 
which  are  the  terms  commonly  used  in  describing  the  operation  of 
ejector-diffuser  systems  (see  Fig.  2).  The  starting  pressure  ratio 
(pex/ Pt^start*  is  the  maximum  exhaust  pressure  ratio  at  which  minimum 
cell  pressure  ratio  is  obtained  for  decreasing  values  of  pex/pt-  The 
operating  pressure  ratio  is  the  maximum  exhaust  pressure  ratio  at  which 
minimum  cell  pressure  ratio  can  be  maintained  for  increasing  values  of 
Pex/Pt-  The  minimum  cell  pressure  ratio  is  obtained  for  a  given  zero- 
secondary-flow,  ejector-diffuser  configuration  when  the  cell  pressure 
becomes  independent  of  exhaust  pressure.  The  starting  pressure  ratio 
is  either  equal  to  or  less  than  the  operating  pressure  ratio  and  is  the 
pressure  ratio  which  is  most  important  in  the  design  of  an  ejector  sys¬ 
tem,  whether  the  pressure  ratio  is  controlled  by  varying  the  exhaust 
pressure  or  by  varying  the  rocket  chamber  pressure.  It  should  be  noted, 
however,  that  minimum  cell  pressure  ratio  is  affected  by  nozzle  total 
pressure  because  of  the  Reynolds  number  effect,  as  previously  reported 
in  Ref.  4.  Since  the  theoretical  calculations  of  (Pc/pt)min  have  been  the 
subject  of  a  previous  investigation  (Ref.  5),  it  will  be  assumed  for  the 
purposes  of  this  report  that  a  method  similar  to  that  discussed  in  Ref.  5 
is  available  for  calculating  this  pressure  ratio. 


SECTION  II 

DEVELOPMENT  OF  THEORY 


The  basic  difference  between  the  method  developed  in  Ref.  2  for  deter¬ 
mining  second-throat  starting  performance  (which  will  be  referred  to 
herein  as  Method  1)  and  the  improved  method  (Method  2)  presented  in  this 
report  is  in  the  method  for  determining  the  conditions  at  station  2  (see 
Fig.  1).  One  of  the  most  important  requirements  for  calculating  the 
ramp  pressure  drag  and  the  resulting  conditions  at  the  second-throat 
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entrance  is  a  knowledge  of  the  driving  nozzle  free-jet  flow  field.  The 
agreement  of  the  experimental  and  theoretical  Mach  number  distribu¬ 
tions  downstream  of  the  nozzle  exit  in  Ref.  5  indicated  that  the  method 
of  characteristics  solution  to  the  general  isentropic  flow  equations  pro¬ 
vides  a  satisfactory  method  for  determining  the  flow  field  approaching 
the  second  throat  for  the  conical  nozzles  investigated  to  date.  As  the 
flow  approaches  the  diffuser  walls,  it  becomes  nonuniform  both  axially 
and  radically.  The  turning  of  such  a  flow  field  parallel  to  the  diffuser 
walls  can  be  treated  approximately  to  determine  the  static  pressure 
distribution  along  the  walls,  as  suggested  in  Ref.  5,  by  considering 
small  increments  in  the  approaching  flow  to  be  uniform  at  some  aver¬ 
age  Mach  number  and  flow  direction.  The  static  pressure  existing 
downstream  of  the  oblique  shock  is  then  assumed  to  be  constant  along 
the  associated  Mach  line  as  shown  in  Fig.  1.  The  Mach  number  preced¬ 
ing  the  ramp  shock  is  determined  from  the  calculated  static  pressures 
on  the  various  Mach  lines  and  the  total  pressure  (which  is  assumed  con¬ 
stant  along  the  streamlines)  downstream  of  the  reflected  boundary  shock. 
A  comparison  of  the  experimental  and  theoretical  pressures  on  the 
second-throat  ramp  in  section  4.  1  indicates  that  good  agreement  results 
with  this  method  until  the  reflected  free-jet  boundary  shock  intersects 
with  the  ramp  shock  (Fig.  1).  The  deviation  at  this  point  is  probably  the 
result  of  neglecting  the  effects  of  boundary  layer  on  the  flow  approaching 
the  ramp  and  the  assumption  that  the  shock  waves  are  not  curved.  An 
empirical  assumption  was  therefore  made  to  correct  for  this  deviation 
by  assuming  that  the  pressure  on  the  second-throat  ramp  is  constant 
downstream  of  the  shock  intersection. 

The  general  procedure  which  is  followed  in  determining  the  condi¬ 
tions  at  the  second-throat  entrance  (station  2)  with  Method  2  can  be 
summarized  as  follows: 

1.  The  minimum  cell  pressure  ratio  is  obtained  using  a 
method  similar  to  that  discussed  in  Ref.  5.  It  is  assumed 
that  cell  pressure  is  not  affected  by  changes  in  diffuser 
geometry  downstream  of  free-jet  impingement. 

2.  The  calculated  cell  pressure  ratio  is  used  with  the  method 
of  characteristic  solution  for  axisymmetric  flow  to  deter¬ 
mine  the  free-jet  flow  field  approaching  the  second  throat. 

3.  The  ramp  pressure  and  the  conditions  entering  the  throat 
are  determined  by  assuming  that  the  flow  is  turned  through 
a  two-dimensional  shock  system  and  that  the  ramp  pressure 
is  constant  after  the  intersection  of  the  reflected  boundary 
shock  and  the  ramp  shock.  The  frictional  forces,  fd  and  fR, 
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the  pressure  drag  force  on  the  second-throat  ramp,  Fr, 
and  the  oblique  shock  location  approaching  the  second- 
throat  entrance  are  then  determined. 

4.  The  results  are  applied  to  the  conservation  equations  to 
obtain  the  conditions  entering  the  second  throat  at  sta¬ 
tion  2,  which  are  then  used  to  determine  the  starting 
pressure  ratio  for  a  given  second-throat  configuration. 


2.1  GENERAL  EQUATIONS 


This  theoretical  analysis  is  a  continuation  of  that  reported  in  Ref.  2 
and  utilizes  the  same  general  equations  for  determining  second-throat 
performance,  which  were  developed  previously  by  applying  the  conserva¬ 
tion  equations  between  stations  1  and  2  and  the  following  assumptions: 

1.  Flow  is  steady,  adiabatic,  and  one -dimensional  at 
stations  2  and  3, 

2.  Gas  is  perfect,  and 


3.  Gas  velocity  in  cell  region  is  zero. 


where 


Fx, 


drag  forces 


(1) 


Fix  —  Fnex  +  Pc  (A(J  Ane ) 


(2) 
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Ld 

fd  =  TTViy  f  C{  pa  M/  dx 
Li 


rd 


Fr  =  2  tt  f  pR  r  d 


rd 


fRx  =  y  cot  @R  f  Cf  Mr  Pr  r  d 

z  rst 

rs  t 

F2x  =  2  77  f  p2  (1  +  y  M22)  r  d  r 


(3) 

(4) 

(5) 

(6) 


By  assuming  the  flow  to  be  one -dimensional  at  station  2,  the  conservation 
"equations  may  be  used  to  write 


1  +  y  M2 


VRTt 


y  —  1  2 

1  +  V- , 


(7) 


Once  the  variables  which  make  up  the  left-hand  side  of  this  equation  have 
been  obtained,  the  Mach  number  may  be  obtained  from  the  graphical 
solution  shown  in  Fig.  3. 


2.2  APPLICATION  OF  THEORY 


The  following  sections  indicate  how  the  theoretical  results  up  to  this 
point  may  be  applied  to  various  second-throat  configurations. 

2.2.1  Application  to  Short  Second  Throats  [  (  L/D)st  <  1  ] 

The  best  method  for  obtaining  the  starting  and  operating  pressure 
ratio  for  short  second-throat  diffusers  is  to  use  the  predetermined  super¬ 
sonic  Mach  number  at  station  2  to  determine  the  maximum  pressure  ratio, 
Pst/P2'  required  to  separate  a  turbulent  boundary  layer  (see  Fig.  1). 

This  value  is  then  used  with  the  continuity  equation 


Pa/Pt 


n/RT 


Pt  Ast  M. 


Vy  ( 


1  +  M, 


(8) 


to  determine  the  overall  starting  pressure  ratio,  pst /  Pt-  This  method 
for  determining  short  second-throat  performance  was  developed  in  Ref.  2 
by  using  the  conservation  equations  and  making  the  additional  assumption 
that  the  frictional  losses  in  the  minimum  area  are  negligible. 

There  is  still  some  question  as  to  which  separation  criteria  should 
be  used.  For  the  semi-empirical  method  (Method  1)  of  determining 
second-throat  performance  developed  in  Ref.  2,  Arens  and  Spiegler's 
separation  pressure  ratio  (Ref.  3)  gave  good  results.  However,  for  the 
theoretical  method  developed  in  this  report  (Method  2),  Mager's  separa¬ 
tion  pressure  ratio  (Ref.  6)  appeared  to  give  better  agreement  with  the 
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diffuser  performance.  The  reason  for  the  disagreement  is  unknown. 

A  comparison  of  the  separation  pressure  ratios  obtained  by  these  two 
methods  is  shown  in  Fig.  4  for  y  =  1.4. 

2.2.2  Application  to  Long  Second  Throats  [  (  L/D  )st  >  5  ] 

The  maximum  diffuser  exit  pressure  at  which  the  long  second- 
throat  diffuser  will  start  and  operate  can  be  determined  by  using  the 
subsonic  solution,  M3,  of  Eq.  (7),  together  with  the  following  addi¬ 
tional  assumptions  for  the  flow  between  stations  2  and  3: 

1.  Frictional  losses  in  the  minimum  area  may  be  estimated 
by  assuming  a  linear  Mach  number  distribution  from 
station  2  to  3  and  a  constant  friction  coefficient,  and 

2.  The  static  pressure  at  station  3  equals  the  exhaust 
pressure,  pex. 

The  starting  and  operating  pressure  can  thus  be  obtained  by  calculating 
the  subsonic  exit  Mach  number,  M3,  and  using  the  previously  calculated 
upstream  Mach  number,  M2,  in  Eq.  (9),  as  developed  in  Ref.  2. 


f'z 


X 


m  yj  R  Tt 


Cf  2  (L/D  )st 
(M3  -  M2 ) 


1  +  y  Mj2 

+  M3Vy(i^M32j 


(9) 


A  graphical  solution  for  M3  is  presented  in  Fig.  3a  for  y  =  1.4, 
(L/D)gt  =  8,  and  Cf  =  0.  003.  The  value  of  M3  is  then  used  in  the  con¬ 
tinuity  equation  to  calculate  the  maximum  starting  and  operating  dif¬ 
fuser  exit  pressure  ratio. 


Pst _ rn  V  R  Tt! _ 

Pt  PtAstM3A/y  (l  +  ^M32j 


(10) 


2.2.3  Applications  to  Intermediate  Length  Second  Throats  [  1  <  (L/D)  <  5  ] 


An  experimental  investigation  of  intermediate  length  second-throat 
performance  indicated  that  the  linear  interpolation  between  the  short  and 
long  second-throat  theories  suggested  in  Ref.  2  gave  the  best  results. 
The  experimental  second-throat  performance  using  conical  nozzles 
(Ref.  2)  and  cylindrical  diffuser  performance  using  annular  nozzles 
(Ref.  7)  both  indicate  that  this  interpolation  can  be  made  between 
(L/D)st  =  0.  5  and  5.  0. 
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An  empirical  mathematical  flow  model  was  also  developed  for  pre¬ 
dicting  the  performance  of  intermediate  length  second  throats.  Although 
this  model  gave  good  results,  it  is  not  a  physically  realistic  flow  model. 
The  development  of  this  mathematical  model  and  its  results  are  included 
in  Appendix  I. 

2.2.4  Application  to  Second  Throats  with  Centerbodies 


The  theoretical  analysis  of  second-throat  diffusers  witn  centerbodies 
uses  the  preceding  theory  with  the  following  additional  assumptions: 

1.  The  pressure  on  the  centerbody  cone  is  calculated  using 
conical  flow  relations,  and 

2.  The  Mach  number  approaching  the  supporting  vanes  is 
subsonic  just  prior  to  unstarting  the  ejector. 


It  should  also  be  noted  that  this  method  has  a  limited  application  in  short, 
second-throat  configurations.  The  centerbody  must  be  positioned  such 
that  the  centerbody  shock  impinges  upstream  of  the  throat  (see  section 
4.2.4). 


STATIONS 


© 


© 


rCB 


F3x  =  F2jc  -  fcB  ~  n  J  PCB  rCB  drcB 


(ID 


where  F2V  is  determined  from  Eq.  (7). 

A 


For  centerbodies  without  supporting  vanes  in  the  throat 

fCB  =  ~  y  Pa  Ma2  n  (Dst  -  Dcb)LcB  ^^) 

and  for  those  using  supporting  vanes  in  the  throat 

fcB  =  “  y  Pa  Ma2  (Dst  -  Dcb)  Lcb  +  2N  (hv)  LVJ  (13) 
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where  Cf  =  0.  003  and  Ma  and  Pa  are  determined  from  the  normal  shock 
relations  for  the  Mach  number,  M2,  which  is  calculated  for  the  condi¬ 
tions  entering  the  throat.  This  Mach  number,  M^,  in  front  of  the  normal 
shock  is  assumed  to  be  an  average  of  the  supersonic  Mach  numbers  down¬ 
stream  of  the  ramp  and  centerbody  shocks  at  station  2.  Equation  (11) 
can  be  further  written  as 

f3x  _  1  +  y  M2 

m  V~R  Tt  Myy^l  +  ^^-M2] 

where  M  is  either  supersonic  or  subsonic,  depending  on  whether  a  short 
or  long  second  throat  is  being  considered.  The  starting  pressure  ratio 
is  then  calculated  using  the  previously  discussed  equations  (Eqs.  (8)  or 
(10))  for  either  a  short  or  long  second  throat. 

Figure  5  shows  that  the  pressure  force  on  the  second-throat  ramp 
as  well  as  the  conditions  entering  the  second  throat  are  affected  at  cer¬ 
tain  locations  of  the  centerbody.  Comparisons  of  theoretical  and  experi¬ 
mental  starting  pressure  ratios  for  the  configurations  investigated  in 
section  4.  2.  4  indicate  that  this  effect  can  be  neglected  providing  devia¬ 
tions  less  than  10  percent  are  not  required. 


SECTION  III 

EXPERIMENTAL  INVESTIGATION 


The  experimental  investigation  of  second-throat  diffuser  perform¬ 
ance  discussed  in  this  report  is  a  continuation  of  the  study  reported  in 
Ref.  2.  Additional  experimental  results  are  presented  in  support  of  the 
following  theoretical  investigations: 

1.  Wall  static  pressure  distribution  on  the  second-throat 
ramp  and  flow  conditions  entering  the  throat, 

2.  Development  of  a  mathematical  flow  model  to  predict 
intermediate  length  second-throat  performance, 

3.  Effects  of  centerbody  position  on  second-throat  per¬ 
formance,  and 

4.  Second-throat  performance  using  annular  and  cluster 
nozzles. 


3.1  EXPERIMENTAL  APPARATUS 

Eight  nozzles  and  eight  second-throat  diffusers  (Table  I)  were  tested 
in  various  combinations  during  this  phase  of  the  investigation  of  the  55 
configurations  listed  in  Table  II. 
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3.1.1  Test  Hardware  Description 


The  nozzles  were  concentrically  located  in  a  cylindrical  diffuser 
with  the  upstream  end  of  the  diffuser  attached  to  a  sealed  plenum.  A 
typical  test  configuration  is  shown  in  Fig.  6.  The  nozzles  were 
mounted  on  a  movable  section  of  inlet  supply  pipe  which  permitted  the 
nozzle  to  be  translated  approximately  9.  0  in.  along  the  horizontal 
centerline  of  the  diffuser  system.  The  design  of  the  "O"  ring  seals 
in  the  telescoping  sections  permitted  the  nozzle  to  be  positioned  during 
a  test  without  leakage  into  the  cell  region.  The  position  of  the  nozzle 
with  respect  to  the  second  throat  was  indicated  by  a  calibrated  counter, 
which  registered  the  rotations  of  the  actuating  mechanism. 

The  nozzle  configurations  which  were  tested  are  illustrated  in 
Fig.  7,  and  the  diffuser  configurations  are  shown  in  Fig.  8.  The  code 
designation  for  the  various  combinations  of  these  configurations  is 
included  in  Table  I.  A  typical  ejector  configuration  designation  would 
be  4-5d-CB,  which  indicates  an  18-deg  conical  nozzle  with  an  area  ratio 
of  25.  0  exhausting  into  a  10.  19-in.  -diam  cylindrical  diffuser  with  a 
second-throat  diffuser  having  a  throat  length  to  diameter  ratio, 

(L/D)st  =  8.  1,  an  effective  contraction  ratio,  (Ast"AcB)/ Ad  =  0.417, 
and  an  18-deg  centerbody  located  in  the  throat. 

3.1.2  Instrumentation 

The  parameters  measured  during  this  investigation  were:  cell  pres¬ 
sure,  pc;  static  pressures  along  the  second  throat  and  cylindrical  diffuser 
walls,  pw;  static  pressures  on  the  second-throat  ramp,  pjj;  nozzle  total 
pressure,  p-^;  nozzle  total  temperature,  T^;  and  second-throat  diffuser 
exit  static  and  total  pressures,  pgt  and  ptr.  Table  III  contains  the  range 
of  the  measured  parameters  and  the  type  of  measuring  instrument  used 
for  each. 


3.2  EXPERIMENTAL  PROCEDURE 

Prior  to  each  test,  the  nozzle,  test  cell,  and  instrumentation  lines 
were  pressure  checked  to  minimize  the  possibility  of  leakage.  A  vacuum 
check  was  also  made  prior  to  each  test  to  reduce  the  possibility  of  instru¬ 
mentation  leakage.  Inlet  air  was  supplied  from  the  RTF  compressors  at 
a  pressure,  p^,  as  high  as  46  psia  and  at  a  temperature  of  approximately 
80°F.  The  ejectors  exhausted  into  the  RTF  exhaust  machines,  which 
provided  pressures  as  low  as  7  mm  HgA.  An  electrically  operated  throt¬ 
tling  valve  in  the  exhaust  ducting  was  used  to  control  the  exhaust  pressure. 
The  inlet  supply  pressure  was  manually  controlled  by  a  gate -type  valve. 
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The  maximum  second-throat  diffuser  exit  pressure,  (Pst)start’  at 
which  the  ejector  became  started  was  obtained  for  each  configuration 
at  a  given  second-throat  and/ or  centerbody  position  and  total  pressure, 
p-j.,  by  decreasing  the  exhaust  pressure  until  the  cell  pressure,  pc, 
reached  a  minimum  value.  The  exhaust  pressure  was  then  increased 
until  the  ejector  again  became  unstarted  (where  pc  started  to  increase) 
to  determine  the  maximum  operating  diffuser  exit  pressure.  This  pro¬ 
cedure  was  repeated  at  various  nozzle  positions  to  determine  the  effect 
of  second-throat  location  on  ejector  second-throat  diffuser  performance. 
The  wall  pressures  on  the  second-throat  diffuser  were  recorded  when 
the  cell  pressure  was  at  the  minimum  value.  During  the  centerbody  in¬ 
vestigation,  the  above  procedures  were  repeated  at  various  nozzle  posi¬ 
tions  and  centerbody  positions  to  determine  the  effects  of  centerbody 
position  on  second-throat  performance.  The  second-throat  exit  total 
pressure,  p-tr,  surveys  made  during  the  intermediate  length  second- 
throat  flow-model  investigation  were  recorded  just  prior  to  unstarting 
the  ejector. 


SECTION  IV 

DISCUSSION  OF  RESULTS 


During  the  previous  second-throat  diffuser  investigation  in  Ref.  2, 
the  empirical  method  for  predicting  the  pressure  force  on  the  second- 
throat  ramp  gave  good  results  for  conical  nozzles  with  uniform  exit 
Mach  numbers.  However,  the  results  shown  in  Table  II  for  annular 
nozzles  indicate  that  this  method  is  not  adequate  for  predicting  the 
starting  pressure  ratio  of  ejector  systems  which  have  highly  nonuni¬ 
form  Mach  number  distribution  at  the  nozzle  exit.  Typical  nonuniform 
nozzle  exit  Mach  number  distributions  are  shown  in  Fig.  9.  Even  the 
flow  from  nozzles  with  uniform  exit  flow  becomes  slightly  nonuniform 
as  it  approaches  the  second  throat  (Fig.  10).  A  more  accurate  proce¬ 
dure  for  predicting  the  pressures  on  the  second-throat  ramp  and  the 
conditions  approaching  the  throat  entrance  (as  discussed  in  section  II) 
was  therefore  developed.  Comparisons  of  experimental  and  theoretical 
results  show  that  this  procedure  is  a  prerequisite  to  the  theoretical 
calculation  of  the  starting  pressure  ratio  for  second-throat  diffusers 
whose  approaching  flow  field  is  highly  nonuniform. 


4.1  RAMP  PRESSURE  FORCE 

Figures  11a  through  f  show  a  comparison  of  the  experimental  and 
theoretical  ramp  pressure  distribution.  It  will  be  noted  that  in  each 
case,  the  theoretical  pressure  distribution  is  less  than  the  experimental 
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value  after  the  intersection  of  the  impingement  shock  and  the  ramp  shock. 
There  are  several  things  which  may  account  for  this.  One  is  the  presence 
of  boundary  layer,  which  was  neglected,  and  another  is  the  assumption 
that  the  impingement  shock  originates  where  the  inviscid  jet  boundary 
strikes  the  cylindrical  diffuser  rather  than  at  the  impingement  of  the  stag¬ 
nating  streamline  as  shown  by  the  following  sketch.  However,  this  devia¬ 
tion  is  more  than  likely  the  result  of  the  two-dimensional  flow 


assumption,  which  becomes  inadequate  as  the  relations  across  the  oblique 
shocks  are  analyzed  at  greater  distances  downstream  of  the  free-jet 
impingement.  The  best  results  for  predicting  the  pressure  force  acting 
on  the  ramp  were  obtained  by  integrating  the  pressure  forces  along  the 
ramp  and  assuming  the  pressure  after  the  shock  intersection  to  be  con¬ 
stant.  Comparisons  of  the  experimental  and  theoretical  ramp  pressure 
forces  using  this  procedure  showed  an  agreement  within  approximately 
five  percent  for  the  conical  nozzle  configurations  investigated.  This 
results  in  an  approximate  1.  7 -percent  deviation  in  the  overall  starting 
pressure  ratio. 

A  theoretical  Method  2  calculation  of  the  ramp  pressures  using 
annular  and  cluster  nozzles  was  not  made  because  of  the  unavailability 
of  a  flow  field.  However,  the  experimental  ramp  pressures  are  shown 
in  Figs.  12a  through  c. 
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Centerbody  position  will  affect  the  pressures  on  the  ramp  as  is 
shown  in  Fig.  5;  however,  a  theoretical  method  for  calculating  this 
effect  was  not  included  in  this  investigation.  Neglecting  the  effect  of 
centerbody  position  in  the  theoretical  calculations  is  reflected  in  the 
percent  deviation  shown  in  Table  II.  A  comparison  of  the  force  on  the 
ramp  for  a  centerbody  located  at  position  (T)  (Fig.  5)  shows  that  the 
force  on  the  ramp  is  increased  approximately  two  and  one-half  times 
that  observed  when  the  centerbody  is  retracted  to  position  (5)  (Fig.  5). 
This  would  result  in  approximately  a  20-percent  deviation  in  the  over¬ 
all  theoretical  starting  pressure  ratio  and  would  indicate  that  the  effect 
of  centerbody  position  should  be  considered  for  an  accurate  calculation. 
However,  later  comparisons  of  the  experimental  and  theoretical  start¬ 
ing  pressure  ratios  (neglecting  the  effect  of  the  centerbody  on  ramp 
pressure)  show  that  the  actual  deviation  was  much  less. 


4.2  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  SECOND-THROAT 
PERFORMANCE 

Method  1,  as  presented  in  Ref.  2,  resulted  in  a  satisfactory  proce¬ 
dure  for  calculating  the  starting  pressure  ratio  for  second-throat 
ejector-diffuser  systems  with  uniform  nozzle  exit  flow.  Although 
Method  2  gives  as  good  or  better  results  than  Method  1  for  these  same 
systems  (see  Table  II),  Method  2  was  developed  for  the  primary  pur¬ 
pose  of  arriving  at  a  more  general  procedure  for  calculating  the  per¬ 
formance  of  systems  whose  nozzle  exit  flow  is  nonuniform.  Unfortunately, 
the  nonuniform  flow  fields  for  the  annular  nozzles  tested  were  not  avail¬ 
able  to  make  a  comparison  of  the  improved  theory  and  the  experimental 
data.  However,  the  good  agreement  of  the  experimental  data  with  the 
results  from  Method  2  for  the  conical  and  contour  nozzles,  whose  flow 
fields  were  available,  should  indicate  that  Method  2  will  be  satisfactory 
once  the  free-jet  flow  field  is  determined.  This  will  be  verified  pending 
further  development  of  a  method  to  calculate  the  annular  nozzle  flow 
fields. 

The  theoretical  starting  pressure  ratio  for  the  annular  nozzles 
listed  in  Table  II  was  calculated  by  using  Method  1  and  assuming  an 
average  nozzle  exit  Mach  number  (see  Fig.  9).  However,  these  results 
are  not  discussed  in  detail  since  Method  1  does  not  provide  an  accurate 
method  for  calculating  the  ramp  force  and  the  conditions  approaching 
the  throat  for  these  nozzles. 

4.2.1  Short  Second  Throats  [  (  L/D  )st  <  1  1 

Comparison  of  the  experimental  and  theoretical  starting  pressure 
ratios  for  six  short  second-throat  configurations  listed  in  Table  II 
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resulted  in  a  maximum  deviation  from  the  experimental  results  of  2.  9  per¬ 
cent  using  Method  2  (see  Fig.  13).  The  new  theoretical  method  predicted 
the  starting  performance  more  accurately  for  all  short  second-throat  con¬ 
figurations  investigated. 

4.2.2  Long  Second  Throots  [  (  L/D  )st  >  5  ] 

The  theoretical  starting  and  operating  pressure  ratio  was  determined 
for  two  long  second-throat  configurations ,  (3-5d  and  4-5d)  using  Method  2. 
Comparisons  with  the  experimental  results  showed  a  maximum  deviation 
of  6.  8  percent  for  conical  nozzles  and  2.  2  percent  for  the  one  contour 
nozzle  investigated  using  Method  2.  Method  1  resulted  in  a  maximum 
deviation  of  9.  7  percent  for  conical  nozzles  and  11.  7  percent  for  the 
contour  nozzle.  Method  2  gives  better  results  for  the  contour  nozzle 
because  it  takes  into  account  the  nonuniform  flow  condition  approaching 
the  second  throat  shown  in  Fig.  10. 

4.2.3  Intermediote  Length  Second  Throots  [  1  <  (  L/D  )st  <  5  ] 

The  starting  pressure  ratio  of  intermediate  length  second  throats  can 
be  determined  within  approximately  10  percent  by  taking  a  linear  inter¬ 
polation  between  the  short  throat  theoretical  results  at  (L/D)st  =  0.  5  and 
the  long  second-throat  theoretical  results  at  (L/D)^  =  5.0  (see  Fig.  14). 

A  more  theoretical  approach  to  determining  the  performance  of  inter¬ 
mediate  length  second  throats  was  attempted,  but  a  satisfactory  proce¬ 
dure  did  not  develop.  The  mathematical  flow  model  which  was  developed 
gave  accuracies  equal  to  the  linear  interpolation;  however,  the  flow  model 
is  questionable.  Comparisons  of  the  experimental  data  with  the  theoretical 
results  for  this  mathematical  flow  model  are  given  in  Appendix  I. 

4.2.4  Centerbody  Second  Throots 

The  performance  of  two  centerbody  second-throat  configurations 
(Fig.  8)  was  investigated  -  one  in  a  short  second  throat  without  supporting 
vanes  in  the  minimum  area  section  and  one  in  a  long  second  throat  with 
supporting  vanes  in  the  throat.  A  comparison  of  the  theoretical  and 
experimental  starting  pressure  ratios  for  the  short  second-throat  con¬ 
figurations  investigated  shows  that  a  maximum  deviation  of  10.  9  percent 
was  obtained  until  the  shock  generated  from  the  centerbody  impinged  on 
the  throat,  as  shown  by  point  (a)  in  Fig.  15a.  The  large  deviation  which 
exists  downstream  of  point  (a)  is  believed  to  be  the  result  of  the  center- 
body  shock  on  the  second-throat  separation  characteristics.  The  theo¬ 
retical  results  upstream  of  point  (a)  were  obtained  using  the  procedures 
outlined  in  Method  2.  However,  as  was  mentioned  in  section  4.  1,  the 
effect  of  the  centerbody  on  second-throat  ramp  pressure  drag  was  not 
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considered.  Figure  15a  further  indicates  that  the  centerbody  configura¬ 
tions  tested  gave  the  best  performance  (maximum  starting  pressure 
ratio)  when  the  centerbody  was  positioned  such  that  the  shock  was 
impinging  upon  the  ramp.  This  is  not  to  be  considered  as  being  generally 
true  since  it  depends  upon  the  geometry  of  the  second  throat  (contraction, 
throat  length,  etc.  ).  When  the  second-throat  contraction  is  near  its 
limit,  the  forces  on  the  ramp  should  be  kept  to  a  minimum. 

Figure  15b  shows  that  for  configuration  3-5d-CB  the  addition  of  a 
centerbody  required  a  decrease  in  pst / Pt  to  start  the  ejector  system. 

An  analysis  of  the  losses  within  the  system  indicated  that  this  was  due 
to  the  large  frictional  losses  caused  by  the  supporting  vanes  in  the  mini¬ 
mum  area  section.  A  comparison  of  the  theoretical  and  experimental 
starting  pressure  ratios  showed  a  maximum  deviation  of  10.  5  percent. 
Figure  15b  shows  that  there  is  a  centerbody  position  where  the  ejector 
would  not  restart.  However,  the  theoretical  calculations  of  the  starting 
performance  for  the  centerbody  in  that  position  indicated  that  the  start¬ 
ing  characteristics  of  the  ejector  were  near  the  starting  limit. 

4.2.5  Cluster  Nozzle  Configurations 


The  second-throat  performance  with  annular  nozzles,  which  has  been 
discussed  in  the  preceding  section,  clearly  indicates  that  a  procedure 
similar  to  the  proposed  method  (Method  2)  is  necessary  to  predict  accu¬ 
rately  the  starting  pressure  ratio.  Cluster  nozzle  configurations  are 
similar  in  that  they,  too,  have  nonuniform  flow  conditions  at  the  nozzle 
exit.  An  attempt  to  determine  theoretically  by  Method  1  the  flow  condi¬ 
tions  approaching  the  second  throat  for  a  four-nozzle  cluster  configura¬ 
tion  shown  in  Fig.  7  was  made  by  using  the  experimental  cell  pressure 
ratio,  an  experimental  average  ramp  pressure,  and  assuming  the  diffuser 
diameter  to  be  =  4.  301  in.  Table  II  indicates  that  a  maximum  deviation 
of  29  percent  results.  This  is  the  result  of  using  an  empirical  procedure 
for  calculating  the  conditions  entering  the  second  throat.  Unfortunately, 
an  adequate  flow  field  was  not  available  to  calculate  theoretically  the 
performance  using  Method  2. 


SECTION  V 

APPLICATION  TO  DESIGN  OF  SECOND-THROAT  DIFFUSER  SYSTEMS 


The  application  of  the  basic  theory  to  given  second-throat  diffuser 
designs  has  been  the  primary  objective  in  the  preceding  sections.  How¬ 
ever,  the  person  who  has  the  problem  of  designing  an  ejector -diffuser 
system  to  successfully  test  a  given  rocket  motor  at  a  specified  altitude 
condition  must  approach  the  theory  from  a  different  aspect.  It  is. 
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therefore,  desirable  to  specify  the  procedures  which  should  be  followed 
to  use  effectively  the  theories  in  the  design  of  a  second-throat  diffuser. 
The  procedures  presented  here  are  general,  such  that  they  can  be  applied 
to  any  of  the  configurations  discussed. 

It  is  first  assumed  that  the  following  information  is  given: 


1. 

Rocket  motor  nozzle  geometry, 

2. 

Ratio  of  specific  heats. 

3. 

Altitude  pressure  desired, 

4. 

Rocket  chamber  pressure. 

5. 

Available  minimum  diffuser  exhaust 

pressure,  and 

6. 

Available  maximum  diffuser  system 

length  and  diameter. 

There  are  several  ways  one  may  approach  the  problem  depending 
upon  the  test  objective.  However,  for  the  purpose  of  discussion,  it  is 
assumed  that  the  main  purpose  of  the  test  is  to  ignite  the  rocket  motor 
and  record  the  steady-state  thrust  at  a  given  altitude. 

The  initial  altitude  condition  prior  to  firing  is  obtained  using  a 
steam  ejector  downstream  of  the  nozzle.  The  details  of  the  design  of 
the  steam  ejector  are  not  discussed  because  they  are  similar  to  the 
design  of  a  hot  rocket  exhaust  gas  ejector  system.  It  should  be  noted, 
however,  that  the  objective  requires  that  the  steam  ejector  obtain  a  cell 
pressure  equal  to  that  which  the  hot  rocket  exhaust  gas  ejector  system 
will  obtain.  If  a  centerbody-type  steam  ejector  is  used,  it  may  be 
desirable  to  locate  the  centerbody  as  close  to  the  nozzle  exit  as  possible 
to  minimize  the  recirculation  of  hot  exhaust  gases  and  variations  in  cell 
pressure  during  the  tailoff  transient  (during  rocket  shutdown).  However, 
the  influence  of  centerbody  position  on  second-throat  performance  (as 
well  as  centerbody  heat  transfer  problems)  must  be  taken  into  account. 
The  configuration  must  be  analyzed  closely  to  ensure  that  the  position 
of  the  centerbody  will  not  interfere  with  starting  of  the  ejector. 

The  first  parameter  to  determine  is  the  capture  diffuser  diameter 
required  to  obtain  the  test  altitude.  This  can  be  found  by  an  iteration 
procedure  using  the  theory  discussed  in  Ref.  5,  which  has  been  pro¬ 
grammed  on  a  7074  computer.  The  flow  field  downstream  of  the  nozzle 
and  the  free -jet  impingement  point  of  the  stagnating  streamline  on  the 
capture  duct  can  then  be  determined  by  using  the  techniques  discussed 
in  Ref.  5.  The  second-throat  ramp  should  be  located  just  downstream 
of  the  free-jet  impingement  point  on  the  cylindrical  diffuser  (see  Ref.  2). 
The  optimum  second-throat  ramp  angle,  which  should  be  selected,  will 
depend  on  the  approaching  flow  conditions;  however,  Ref.  1  indicated 
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that  the  starting  pressure  ratio  was  rather  insensitive  to  ramp  angle 
for  the  conical  nozzles  investigated.  (Ramp  angle  was  varied  from 
6  to  18  deg. ) 

The  next  step  would  be  to  determine  the  pressure  on  the  second- 
throat  ramp.  To  accomplish  this,  the  steps  outlined  in  the  beginning 
of  section  II  are  required.  The  maximum  second-throat  contraction 
can,  therefore,  be  approximated  by  assuming  M2  =  1.0  and  solving 
for  F2x  in  Eq.  (7).  The  second-throat  contraction  can  then  be  calcu¬ 
lated  using  Eq.  (1)  and  neglecting  the  frictional  forces  for  the  first 
approximation.  If  a  centerbody  is  employed  in  the  second  throat,  then 
several  iterations  are  required  to  determine  the  values  of  Ma  and  Pa 
in  Eq.  (12)  or  (13)  since  the  centerbody  frictional  forces  may  become 
appreciable.  Also,  the  design  contraction  ratio,  Agt/A^,  should  have 
at  least  a  10 -percent  margin  of  safety  greater  than  the  limiting  contrac¬ 
tion  ratio  to  ensure  that  the  ejector  will  start.  Once  the  design  contrac¬ 
tion  is  determined,  the  length  of  the  minimum  area  section  should  be 
made  as  long  as  possible  (up  to  approximately  five  throat  diameters). 
This  analysis  also  assumes  no  subsonic  diffuser  is  used.  The  maximum 
starting  pressure  can  then  be  calculated  using  one  of  the  applicable  equa¬ 
tions  given  in  the  preceding  sections.  Should  the  calculated  starting 
pressure  ratio  be  less  than  the  available  starting  pressure  ratio,  it  may 
be  necessary  to  reduce  the  capture  diffuser  diameter  and  compromise 
the  altitude  requirements  to  obtain  a  system  which  will  operate  at  the 
available  pressures. 


SECTION  VI 
CONCLUSIONS 


Based  on  the  results  of  this  investigation,  it  can  be  concluded  that 
an  accurate  theoretical  method  has  been  developed  for  analyzing  the 
starting  characteristics  of  a  given  second-throat,  ejector-diffuser  sys¬ 
tem  subject  to  the  following  limitations: 

1.  The  location  of  the  free-jet  impingement  of  the  stagnating 
streamline  must  be  determined, 

2.  The  flow  field  approaching  the  second  throat  must  be 
known,  and 

3.  For  short  second-throat  diffusers  with  centerbodies, 
the  shock  generated  by  the  centerbody  must  impinge 
upstream  of  the  minimum  area  section. 
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Further  work  which  should  be  accomplished  to  understand  more 
fully  some  of  the  observed  experimental  trends  includes 

1.  The  development  of  a  computer  program  for  determining 
the  flow  field  approaching  and  entering  the  minimum  area 
section  of  a  second  throat  which  could  be  applied  to  noz¬ 
zles  whose  exit  flow  field  is  nonuniform,  and 

2.  A  study  to  determine  the  internal  flow  separation  phenomena 
which  take  place  in  intermediate  length  second  throats  and 
in  short  second  throats  with  a  centerbody  positioned  such 
that  the  generated  shock  impinges  in  the  throat. 
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Fig.  1  Typical  Flow  Model  for  Second-Throat  Operation 
Just  Prior  to  Unstart  Using  a  Conical  Nozzle 


Pex/pt 

Fig.  2  Typical  Ejector  System  Starting  Phenomena  for  Constant 
Nozzle  Total  Pressure 
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Fig.  6  Typical  Ejector-Diffuser  Configuration 


Expansion-Deflection 
(E-D)  Nozzle 


U-  2.36-J 


d.  Forced-Deflection  (F-D)  Nozzle 


AE  DC-TR-65- 1 24 


b.  Centerbody  Test  Configuration  (Dd  =  6.01  in.) 


Fig.  8  Second-Throat  Diffuser  Configurations  Investigated 
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NOZZLE  UP 
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b.  Annular  Spike  Nozzle  (S2-P1) 

Fig.  9  Nozzle  Exit  Mach  Number  Distribution  for  Annular  Nozzles 
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Fig.  10  Mach  Number  Distribution  at  Second-Throat  Ramp  Lip  Station 
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a.  Configuration  3-5d,  (L/D)j  -  0.80 

Fig.  11  Experimental  and  Theoretical  Static  Pressure  Distributions  on  Second-Throat  Ramp  (for  Minimum  Pex/pt) 
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Fig.  11  Continued 
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Configuration  3-6,  (L/D)j 
Fig.  11  Continued 


Fig.  11  Continued 


b.  Configuration  S2P1-5a 


c.  Configuration  9-6 


Fig.  12  Concluded 
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Fig.  13  Comparison  of  Experimental  and  Theoretical  Second-Throat 
Starting  Pressure  Ratio  (Configuration  4-6) 
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Fig.  14  Comparison  of  Experimental  and  Theoretical  Intermediate  Length  Second-Throat  Starting  Pressure  Ratio 
Using  Linear  Interpolation 


AEDC-TR-65-124 


TABLE  I 

DESCRIPTION  OF  NOZZLES  AND  SECOND  THROATS 


o.  Nozzl 

e  Dimensions 

Config. 

Code 

Ane/A* 

A*, 

in.  2 

e 

deg 

Dne. 

in. 

Type  Nozzle 

Ref. 

Fig.  No. 

3+ 

10.  85 

1.  251 

18 

4.  16 

Conical 

7a 

4+ 

25.  00 

0.  542 

18 

4.  16 

Conical 

7a 

5+ 

23.  68 

0.  636 

0 

4.  38 

Contour 

7b 

8+ 

18.  57 

0.  754 

18 

4.  21 

Conical 

7a 

9+ 

10.  7++ 

0.  305 

7.  5++ 

1.  00++ 

Four-Nozzle  Cluster 

7f 

E-D 

20.  0 

0.  937 

8 

4.  196 

Annular 

7c 

F-D 

21.  1 

0.  654 

8 

4.  196 

Annular 

7d 

S2-P1 

17.  5 

0.  564 

18-5 

3.  546 

Annular 

7e 

b. 

Second- Throot  Diffuser  Dimensions 

Config. 

Ast/ Ad 

0R. 

Ad, 

(L/D)st 

Ast  -  ACB 

®CB 

_ 

Ref. 

Code 

deg 

in.  2 

Ad 

Fig.  No. 

5a 

0.  438 

12 

81.  55 

0.  55 

— 

8a 

5e 

0.438 

12 

81.  55 

1.5 

— 

— 

8a 

5b 

0.  438 

12 

81.  55 

3.  0 

— 

— 

8a 

5c 

0.  438 

12 

81.  55 

5.0 

— 

— 

8a 

5d 

0.  438 

12 

81.  55 

8.  1 

— 

— 

8a 

5d-CB 

0.  438 

12 

81.  55 

8.  1 

0.  417 

18 

8c 

6 

0.  576 

6 

28.  31 

0.  21 

— 

— 

8a 

6-CB 

0.  576 

6 

28.  31 

0.  21 

0.  482 

18 

8b 

+  Same  code  as  Ref.  2 
++  Description  for  each  nozzle 
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TABLE  II 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

a.  Long  Second-Throat  Performance 


Config. 

Code 

Ane 

A* 

Od> 

in. 

Ast 

Ad 

(L/D)st 

/l\ 

Wd 

psia 

_ 

Exper. 

(Pc/Pt) 

Exper. 
<Pst  /  Pt> 
Start 

Exper. 
<Pst/  Pt) 
Oper. 

Method  1 

Method  2 

Theory 
(Pst/ Pt) 
Start 

Percent 

Deviation 

Theory 

(Pst/Pt) 

Start 

Percent 

Deviation 

3-5d 

10 

85 

10.  19 

0.438 

8. 

1 

0.  74 

45 

3 

0. 00056+ 

0. 0230 ■ 

0. 0280+ 

0. 0252 

-  9.  7 

0.  0299 

+  6.  8 

3-5d 

10 

85 

0.  80 

45 

3 

0.  00056+ 

0.  0285+ 

0. 0285+ 

0. 0288 

+  1.  1 

0.  0297 

+4.2 

4-5d 

25.  0 

0.  80 

44 

9 

0.  00293+ 

0. 0136+ 

0. 0136+ 

0. 0128 

-  5.  9 

0.0138 

+  1.5 

ED- 

5d 

20 

0 

0.  67 

45 

0 

0.000223 

0.  0173 

0.  0173 

0. 0144 

-16.  8 

— 

— 

0.  80 

0. 000224 

0.0197 

0.0197 

0. 0156 

-20.  8 

— 

— 

0.  92 

0.000224 

0.0199 

0. 0199 

0.0162 

-18.  6 

— 

- 

ED- 

3d 

20 

0 

1.  17 

45 

0 

0. 000224 

0. 0189 

0.0189 

0.0175 

-7.  4 

— 

- 

3d 

23 

68 

0.  80 

44 

8 

0.  000155 

0.  0169 

0.  0171 

0.  0181 

+  7.  1 

- 

l  ^ 

0.  86 

45 

0 

0. 000155 

0.0167 

0.0169 

0. 0181 

+  8.4 

— 

_ 

5- 

3d 

23 

68 

10.  19 

0.438 

8 

1 

0.92 

45 

0 

0.000151 

0.0162 

0.  0167 

0.0181 

+  11.  7 

0.0166 

+2.2 

b. 

Short  Second-Throat  Performance 

3-5a 

10 

85 

10.  19 

0. 

138 

0 

55 

0.  74 

44.  9 

0. 000551+ 

0.0169+ 

0. 0175+ 

0.0165 

-2.  4 

0.  0166 

-1.8 

3-5a 

10.  85 

0.  80 

44 

.9 

0.000551+ 

0. 0170+ 

0.  0177+ 

0.0150 

-13.  3 

0.0168 

-1.2 

4- 

5a 

25.0 

0.  80 

44.5 

0.000299+ 

0.  00648+ 

0.00665+ 

0.  00636 

-1.8 

0.00656 

+  1.2 

FD-5a 

21 

1 

0.  68 

44 

.3 

0. 00016 

0. 0052 

0. 0068 

— 

- 

- 

— 

FD-5a 

21 

1 

0.  92 

44 

4 

0. 00016 

0. 0053 

0.0059 

— 

- 

- 

— 

!S2P1 

-5a 

17 

5 

0.  63 

45 

.4 

0. 000280 

0.00531 

0.  00543 

0. 00580 

+  9.  0 

— 

- 

1 

0.  74 

45 

1 

0. 000276 

0.00597 

0.  00598 

0. 00644 

+  7.  9 

- 

_ 

S2P1 

-5a 

17 

5 

10 

19 

0.438 

0 

55 

1.41 

45 

4 

0.000240 

0.00356 

0.  00367 

0.  00541 

+52.  0 

- 

3- 

S 

10 

85 

6 

01 

0. 

)77 

0 

21 

0.  60 

44.  9 

0.00204 

0. 0263 

0.0298 

0. 0270 

+2.  6 

- 

- 

3- 

S 

10 

.85 

0.  70 

44 

8 

0.00188 

0.  0275 

0.  0275 

0.  0268 

-2.  6 

0. 0283 

+2.  9 

4- 

S 

25 

.0 

0.  60 

45 

7 

0. 000888 

0.0119 

0.  0119 

0. 01062 

-12.0 

0.  0120 

+0.  8 

4- 

S 

25 

0 

0.  80 

45 

5 

0. 000891 

0. 0107 

0.  0107 

0. 01025 

-4.  3 

0. 01092 

+2.  0 

Lt 

S 

10 

7 

6 

01 

0. 

377 

0 

21 

1.  02 

43 

.3 

0.000392 

0. 00974 

0.01095 

0.  00692 

-29.0 

J 

— 

+Data  taken  from  Ref.  2 
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TABLE  li  (Continued) 

c.  Intermediote  Length  Second- Throot  Performonce  (Dj  =  10.19,  A^/A^  =  0.438) 


Config. 

Code 

(L/D)st 

<L/D)d 

Pt- 

psia 

Exper. 

(pc/pt) 

Exper. 

(Pst/r>t) 

Start 

Exper. 
(Pst/  Pt) 
Oper. 

Exper. 

<Pex/  Pt> 
Start 

Exper. 

<P  ex/  Pt) 
Oper. 

Method 

1 

Method 

2 

Theory 

<P  ex/  Pt) 
Start 

Percent 

Deviation 

Theory 
(Pex/  Pt> 
Start 

Percent 

Deviation 

4-5 

e 

1.  50 

1.  18 

44.  12 

— 

0. 00510 

0. 00596 

0.00633 

0.  00664 

0. 00579 

-9.  3 

— 

— 

1.06 

44.  12 

— 

0.00568 

0.00580 

0.00723 

0.00735 

0.00625 

-13.5 

— 

— 

0.  94 

44.  12 

— 

0. 00658 

0. 00718 

0.00788 

0.  00815 

0.  00713 

-9.  5 

— 

— 

4-5e 

i. 

0 

0.82 

44.02 

0.00830 

0. 00850 

0.00857 

0.00892 

0. 00826 

-3.  6 

0.00950 

+10.40 

4-5 

b 

3. 

)0 

1.  18 

44.07 

0.000280 

0. 00663 

0.00714 

0.00878 

0.  00900 

0. 00975 

+  11.0 

— 

— 

0.  94 

44.33 

0.  000274 

0.  01167 

0. 01177 

0. 01090 

0.01130 

0. 01070 

-1.  8 

- 

- 

4-5 

b 

3. 

)0 

0.  82 

44.23 

- 

0.  00933 

0.00991 

0.01240 

0.  01250 

0. 01146 

-7.3 

0.  0129 

+4.  00 

4-5c 

5. 

)0 

1.06 

44.  23 

0.  000276 

0.01322 

0.0134 

0.01350 

0.01360 

0. 01411 

+4.  5 

_ 

— 

0.  94 

44.  12 

0.  000268 

— 

0.  01370 

0.01370 

0.01397 

+2.  0 

_ 

— 

DO 

0.  82 

44.43 

0.  000267 

- 

- 

0. 01350 

0.  01350 

0. 01381 

+2.  2 

0.01372 

+  1.  60 

}  F-D-5b 

3. 

30 

1.  18 

42.93 

0.  000207 

— 

— 

0.  01380 

0.  01400 

— 

— 

, 

0.  94 

42.  93 

0.  000207 

_ 

0.  01390 

0.01510 

— 

- 

— 

- 

0.  82 

42.  93 

0. 000198 

- 

— 

0. 01320 

0. 01650 

— 

- 

— 

—  1 

0.  79 

42.  93 

0.  000166 

— 

— 

0. 01320 

0.  01590 

- 

- 

— 

— 

F- 

D-5b 

3. 

30 

0.  70 

42.  93 

0.  000133 

— 

- 

0. 01300 

0. 01660 

- 

- 

- 

S2-Pl-5b 

3.  00 

0.  36 

44.  32 

0. 000292 

— 

0.01000 

0.01040 

- 

— 

| 

S2-Pl-5b 

3.  00 

_ 

0.  23 

_ 

44.42 

0.  000294 

— 

- 

0. 01090 

0.01110 

- 

8- 

5b 

3. 

DO 

1.42 

45.  52 

0. 000360 

— 

— 

0.01280 

0.01350 

— 

— 

1.  18 

45.  42 

0. 000345 

_ 

0.01430 

0.01440 

- 

_ 

- 

0.94 

45.42 

— 

_ 

0.01460 

0. 01500 

— 

— 

- 

— 

0.  82 

45.42 

0. 000361 

0. 01650 

0.  01670 

0.  01575 

-4.  5 

— 

- 

8- 

5b 

3. 

30 

0.  70 

45.37 

0. 000353 

0.01580 

0.01600 

0. 01572 

-0.  5 

- 

— 

00 
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op 


4-5d-CB  25.00 


Ast/ -^d  (L/D)st 


0.438 

1 

;  8.  1 

.0 

0.438 

8.  j 

10 

0.  577 

0.2 

>1 

0.  577 

0.  21 

TABLE  II  (Cone 

luded) 

d.  Centerbody  Perfo 

rmance 

] 

L 

Method  2 

1  Ast  -  ACB  XCB 
Ad  Dd 

Pt, 

;  psia 

i  Exper.  - 
<pc/  pt) 

Exper. 

(pst/pt) 

Exper.  , 
(Pst/Pt)  | 

Theory  1 

(Pct/Pt)  1  Percent 
Deviation 

45.2  — 

45.2  0.00055 

44.  6  — 

!  44.8  0.00030 

44.8  0.00190 

1  45.4  0.00189 

145.0  0.  00184 

'45.  7  0.00090 

I 

45.2  0.00090 

145.6  0.00090 


LIMIT 

Ejec.  will  not  restart 

0. 0228  I  0. 0228 


LIMIT 

(Ejec.  will  not  restart 
I 

0.02280  0 
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TABLE  III 

DESCRIPTION  OF  MEASURING  INSTRUMENTS 


Parameter 

Measured 

Range 

Measured 

Measuring 

Instrument 

Pc 

0.2  to  5  mm  HgA 

McLeod 

(with  Nitrogen  Cold  Trap) 

5  to  50  mm  HgA 

Diaphragm -Activated 
Dial  Gage 

Pex,  Pst,  and 

Ps 

7  to  50  mm  HgA 

Diaphragm -Activated 
Dial  Gage 

1  to  10  psia 

Diaphragm -Activated 
Dial  Gage 

Pt 

1  to  46  psia 

Diaphragm -Activated 
Dial  Gage 

Pw,  PR>  and 

Ptr 

0.  1  to  90  in.  (Oil) 

Manometer  (Silicone 
Oil-sp  gr  =  1.  092  at  80°F) 

Tt 

70  to  100°F 

Copper -Constantan 
Thermocouple 

Ptr 

0.2  to  33  in.  HgA 

Manometer  (Mercury) 
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APPENDIX  I 

INTERMEDIATE  LENGTH  SECOND-THROAT  MATHEMATICAL 
FLOW  MODEL  DEVELOPMENT 


Because  of  the  complexity  of  the  intermediate  length  diffuser  flow 
field,  a  mathematical  flow  model  was  developed  to  aid  in  calculating  the 
starting  pressure  ratio  for  intermediate  length  second-throat  diffusers. 
This  model  is  essentially  an  interpolation  of  the  two  flow  models 
assumed  for  short  and  long  second  throats.  The  following  sketch  shows 
typical  flow  models  which  could  be  assumed  for  long  and  short  second- 
throat  diffusers. 


The  long  second-throat  model  assumes  the  flow  shocks  down,  with  one- 
dimensional  subsonic  flow  existing  at  station  3.  The  short  second- 
throat  model  assumes  a  normal  shock  at  station  2.  Both  flow  models 
assume  separation  of  the  turbulent  boundary  layer  along  the  diffuser 
walls.  These  two  flow  models  were  substantiated  only  partially  experi¬ 
mentally  and  are,  therefore,  offered  as  one  possible  description  of  the 
phenomena  that  actually  take  place. 


INTERMEDIATE  LENGTH  SECOND-THROAT  RAKE  STUDY 

The  total  pressure  profile  was  obtained  at  the  exit  of  the  second- 
throat  diffuser  in  an  attempt  to  learn  more  about  the  flow  structure  at 
station  3  for  various  throat  lengths,  as  shown  in  the  sketch  on  the 
following  page. 
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Typical  Second-Throat  Rake  Installation 

Figure  1-1  shows  the  total  pressure  profiles  which  were  taken  at  the 
exit  of  second-throat  configurations  4-5a,  -5b,  -5c,  and  -5e,  at  various 
nozzle  positions. 

There  are  two  methods  for  calculating  the  stream  Mach  number. 

One  method  uses  the  measured  static  and  rake  total  pressures,  and 
the  other  uses  the  measured  rake  and  nozzle  total  pressures.  Whether 
subsonic  or  supersonic  flow  exists  can  definitely  be  determined  by  the 
first  method.  However,  it  is  impossible  to  determine  from  total  pres¬ 
sure  measurements  alone  if  the  measured  total  pressure  loss  was  due 
to  a  normal  standing  in  the  second  throat  (subsonic  solution)  or  to  a 
probe -induced  shock  (supersonic  solution).  An  examination  of  the  meas¬ 
ured  total  pressure  profiles  shown  in  Fig.  I- la  shows  relative  uniform 
profiles  for  second-throat  lengths  of  (L/D)g^-  =  5  and  3.  It  should  be 
noted  that  the  pressure  profiles  shown  in  Figs.  I- la  and  b  were  made 
just  prior  to  unstarting  the  ejector  system.  An  assumption  was,  there¬ 
fore,  made  that  the  measured  wall  static  pressure  for  these  configura¬ 
tions  was  equal  to  the  stream  static  pressure  at  this  station.  The  calcu¬ 
lated  pressure  ratio,  Pst/ptr>  indicated  that  the  flow  was  definitely 
subsonic.  However,  when  rake  measurements  were  made  for  second- 
throat  lengths  of  (L/D)s^-  =  1.5  and  0.  55,  a  definite  change  in  the  profile 
was  noted,  and  the  stream  static  pressure  could  not  be  assumed  equal 
to  the  wall  static  pressure.  Although  it  cannot  be  determined  definitely 
that  supersonic  or  subsonic  flow  exists  at  station  3  from  the  rake  data, 
it  is  feasible  that  a  subsonic  core  could  exist  for  the  short  throat  con¬ 
figuration. 


FLOW  MODEL  EQUATIONS 


A  mathematical  flow  model  was  developed  based  on  the  experi¬ 
mental  observations  and  the  premise  that  the  subsonic  flow  could  exist 
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at  station  3  for  both  long  and  short  second-throat  diffusers. 


STATIONS  ©  ® 


This  model  assumes  a  uniform  subsonic  core  surrounded  by  a  uni¬ 
form  supersonic  annulus  with  both  regions  at  the  same  static  pressure. 

It  should  be  noted  that  the  various  intermediate  second-throat  lengths 
probably  do  not  have  a  pure  subsonic  core  just  prior  to  starting;  how¬ 
ever,  the  model  merely  assumes  the  subsonic  Mach  number  which 
would  exist  at  station  3  if  the  flow  were  to  shock  down  between  stations 
2  and  3  for  each  intermediate  length. 

By  making  additional  assumptions  of 

1.  Steady  flow, 

2.  Perfect  gas,  and 

3.  Adiabatic  flow, 

the  conservation  equations  may  be  written  as 

F3x _ Ac(l+yMc2)  +  (A.t- Ac)  (l  +  y  Ma* ) _ 

m  VRTt  Ac  McVy  (l  +  Mc2)  +  (A„t- Ac)  M„VV  (l  1 

In  addition,  it  is  assumed  that  the  subsonic  core  is  formed  by  a 
normal  shock  at  station  2,  the  size  of  which  is  determined  by  the  posi¬ 
tion  of  the  ramp  shock  at  this  station.  Both  the  position  of  the  ramp 
shock  at  station  2  and  F3X  can  be  obtained  from  the  previous  theoretical 
calculations  of  the  forces  upstream  of  the  throat  entrance  using  Method  2. 
The  subsonic  core  Mach  number,  Mc,  at  station  2  can  be  found  by  the 
normal  shock  relations  for  the  existing  supersonic  free-stream  center- 
line  Mach  number,  Mq^,  at  station  2.  This  allows  the  solution  of  the 
supersonic  Mach  number  in  the  annulus,  Ma>  station  2.  If  the  flow 
field  approaching  the  second  throat  is  not  available,  the  following  section 
provides  an  alternate  semi-empirical  solution  which  works  reasonably 
well  for  nozzles  with  uniform  exit  conditions. 
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In  order  to  solve  Eq.  (1-1)  at  station  3,  two  additional  assumptions 
are  made. 

1.  A  linear  decrease  of  Ma  from  the  calculated  value 
at  station  2  to  a  value  of  unity  of  five  diameters 
downstream  and 

2.  A  linear  growth  of  Ac  from  the  calculated  value  at 
station  2  to  Ast  at  five  diameters  downstream. 

The  maximum  starting  diffuser  exit  pressure  ratio  for  any  given 
second-throat  length  between  (L/D)s^  =  1  and  5  may  now  be  calculated 
from 

Pst  _  _ m  \J  RTt/pt _ 

Pt  ACJ  M C3Vy  (l  +  Mc32J  +  (Ast-A03)  Ma3  Vy  (l  +  rrMajJ) 


SEMI-EMPIRICAL  METHOD  FOR  DETERMINING  CONDITIONS  AT  THROAT  ENTRANCE 


Both  the  size  of  the  normal  shock  disk  assumed  at  station  2  and  the 
Mach  number  approaching  this  disk  must  be  determined  in  order  to 
solve  the  intermediate  length  second-throat  starting  model.  The  core 
Mach  number  approaching  this  disk  may  be  rigorously  determined  by  a 
straight  forward  application  of  an  irrotational  inviscid  characteristics 
solution.  If  a  conical  nozzle  is  used  as  the  driving  nozzle,  the  center- 
line  Mach  number  distribution  may  rigorously  be  determined  from  a 
simple  closed  form  solution.  This  method  is  based  on  the  fact  that  the 
free -jet  core  flow  is  conical  in  nature  until  the  first  expansion  line  from 
the  nozzle  lip  crosses  the  nozzle  centerline.  For  a  high  area  ratio 
nozzle,  this  point  is  generally  downstream  of  station  2.  The  fact  that 
the  centerline  flow  is  conical  resolves  to  the  following  equation: 


A* 


(1-3) 


where  Ax/A*  is  the  expansion  area  ratio  experienced  by  the  gas  on  the 
nozzle  centerline  and  X  is  the  distance  downstream  of  the  nozzle  exit. 
The  centerline  Mach  number  distribution  may  then  be  determined  using 
isentropic  compressible  flow  tables. 


Location  of  the  ramp  shock  at  station  2  is  considerably  more  diffi¬ 
cult  because  the  nonuniform  flow  field  from  the  nozzle  at  this  station 
causes  the  flow  to  be  rotational  behind  the  shock.  A  semi-empirical 
solution  of  the  shock  location  at  station  2  assumes  the  ramp  shock  is 
two-dimensional  (straight)  and  emanates  half-way  between  the  jet 
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impingement  point  and  the  ramp  leading  edge.  It  is  further  assumed 
that  the  gas  expands  from  the  nozzle  to  a  mean  gas  expansion  angle 
(0m)  and  is  turned  parallel  to  the  second-throat  ramp  through  the  two- 
dimensional  oblique  shock  system.  This  mean  gas  expansion  angle 
(0m)  was  developed  in  Ref.  2  for  estimating  ramp  pressure  area  drag 
and  may  be  calculated  from  the  following  equation: 


—  1 

dm  =  tan 


tan  dst 


+  1 


d-4) 


The  Mach  number  used  with  the  mean  gas  expansion  angle  is  the  isen- 
tropic  value  based  on  the  area  ratio,  A^/A*. 


COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  INTERMEDIATE  LENGTH 
SECOND-THROAT  PERFORMANCE 

The  majority  of  the  investigation  of  a  theoretical  method  for  calcu¬ 
lating  intermediate  length  second-throat  performance  was  accomplished 
using  Method  1,  as  is  evident  from  the  comparisons  of  theoretical  and 
experimental  starting  pressure  ratios  in  Figs.  I-2a  and  b  and  as  tabu¬ 
lated  in  Table  II.  Figure  I-2a  compares  the  results  using  Method  1  and 
the  intermediate  length  second-throat  flow  model  with  the  available 
conical  nozzle  starting  pressure  ratio  data  for  optimum  second-throat 
positions.  If  the  gross  simplification  of  the  flow  model  is  considered, 
good  agreement  was  obtained  between  the  experimental  and  theoretical 
results. 

It  is  anticipated  that  Method  2  will  be  required  for  nozzles  whose 
exit  Mach  number  distributions  are  nonuniform.  Figure  I-2a  shows  a 
comparison  of  both  methods  for  the  Ane/A*  =  25.0  conical  nozzle. 

Both  Method  1  and  2  agreed  with  the  experimental  results  within  approxi¬ 
mately  10  percent.  Table  II  shows  that  the  percentage  varies  with  the 
second-throat  length. 

Figure  I-2b  compares  the  theoretical  and  experimental  second-throat 
starting  pressure  ratios  at  various  second-throat  positions  for  the 
Arne/A.*  =  25.  0  conical  nozzle  using  Method  1.  The  maximum  error 
(12 -percent  deviation)  was  found  at  the  maximum  Xst/Dd.  The  maximum 
error  at  this  point  is  thought  to  be  the  result  of  the  assumptions  that  the 
normal  shock  disk  is  located  at  station  2  and  in  the  method  used  to  calcu¬ 
late  the  size  of  this  disk.  The  assumption  of  uniform  static  pressure  at 
each  station,  also,  is  less  valid  at  minimum  (L/D)s^..  Because  of  these 
discrepancies  in  the  assumptions  at  small  values  of  (L/D)st,  it  is  recom¬ 
mended  that  the  short  second-throat  theory  described  in  section  2.  2.  1  be 
used  for  (L/D)st  less  than  one. 
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(PeAt> 


b.  Effect  of  Second-Throat  Position 


Fig.  1-2  Concluded 
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